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Abstract 

Using the approach of autonomous Hamiltonian dynamical system, we at- 
tempt to estimate the as yet unknown upper limit on the size of the galactic 
halo surrounding galaxies (lenses). The key to determine the size of the halo is 
to determine the maximum radius up to which stable material circular orbits 
are admissible. We shall illustrate the approach by considering a solution of 
the Weyl gravity containing a halo parameter 7. The upper limit for several 
observed lenses are calculated for a typical value of 7 for definiteness, with and 
without the cosmological constant A. These lenses (all having Einstein radius 
Re ~ 10^'^ cm) consistently yield an upper limit i?^|!j'°(~ 4.25 x 10^'' cm) inside 
the de Sitter radius only when 7 is negative, while a positive 7 yields -Rmax'** 
always exceeding the de Sitter radius. 
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The issue of dark matter (halo matter), arising out of reconciling known 
gravitational laws with observed flat rotation curves, is a challenging problem in 
modern astrophysics. Observationally, so far, mass of the dark matter appears to 
increase with distance in galaxies, but in clusters exactly the reverse is true, the 
dark matter distribution actually decreases with distance. Indeed, for certain 
dwarfs (such as DD0154) the rotation curve has been measured to almost 15 
optical length scales indicating that the dark matter surrounding this object 
is extremely spread out (see, for a review, Sahni [1]). The total mass of an 
individual galaxy is still somewhat of an unknown quantity since a turn around 
to the V oc r~^/^ law at large radii has not been convincingly observed. On 
the other hand, dark matter is attractive and localized on scales less than the 
cosmological distances where repulsive dark energy prevails. Therefore, one 
would like to know if there could be any upper limit on the size of the dark 
matter halo. The motivation here is to investigate this question using a certain 
model theory. 

Although classical Einstein's general relativity theory has been nicely con- 
firmed within the weak field regime of solar gravity and binary pulsars, observa- 
tions of flat rotation curves in the galactic halo still lack a universally accepted 
satisfactory explanation. The most widely accepted explanation hypothesizes 
that almost every galaxy hosts a large amount of nonluminous matter, the so 
called gravitational dark matter, consisting of unknown particles not included 
in the particle standard model, forming a halo around the galaxy. One of the 
possibility could be that these particles (WIMPs) encircling the galactic center 
are localized in a thick shell providing the needed gravitational fleld [2]. The 
exact nature of either the dark matter or dark energy is yet far too unknown be- 
yond such speculations. There also exist alternative theories, such as Modified 
Newtonian Dynamics (MOND) [3,4], braneworld model [5], scalar fleld model 
[6] etc that attempt to explain dark matter without hypothesizing them. A 
prominent model theory is Weyl conformal gravity and a particular solution in 
the theory is the Mannheim-Kazanas-de Sitter (MKdS) metric [7] that we shall 
consider here. 

The MKdS solution contains two arbitrary parameters 7 and k (= A/3) that 
are expected to play prominent roles on the galactic halo and cosmological scales 
respectively. While the value of the cosmological constant A = 1.29x10"^^ cm~^ 
is well accepted, there is some ambiguity about the sign and magnitude of 7. 
From the flat rotation curve data, Mannheim and Kazanas flx it to be positive 
and being of the order of the inverse Hubble length, while Pireaux [8] argues for 
I7I ~ 10~^^ cm~^. Edery and Paranjape [9] obtained a negative value from the 
gravitational time delay by galactic clusters while the magnitude is still of the 
order of inverse Hiibble length. 

In the present context, we recall that massive neutral hydrogen atoms are ex- 
ecuting circular motions in the halo around the galactic center. The redshifted 
light from those atoms are measured to determine their tangential velocities 
[10,11]. Therefore it is necessary to consider massive test particle orbits. How- 



2 



ever, due to conformal invariance of the theory, geodesies for massive partielcs 
would in general depend on the conformal factor fl'^{x), but here we assume 
a fixed conformal frame and not considering other conformal variants of the 

metric. 

The purpose of this article is to determine the maximum radius within which 
there can exist stable circular orbits of massive test particles, beyond which the 
orbits become unstable. The criterion of stability has been originally suggested 
by Edery and Paranjape [9] because it provides a way to the determination of a 
natural length scale or region of influence of localized sources in the cosmological 
setting. The strategy we adopt here is to frame the geodesic equation in the 
MKdS solution as a Hamiltonian system, and based on it, analyze the stability, 
which seems to favor a negative 7. 

The MKdS metric is given by [7,9] (vacuum speed of light, cq = 1): 



(1) 



dT^ = B{r)dt'^-^dr'^-r^{de^+sm^0dip^), B{r) = 

where k and 7 are constants. 

Using u = 1/r, we get the following path equation for a test particle of mass 
mo on the equatorial plane 6 = Tr/2: 

d'^u 2 7 ^ 1 / 2k\ 

where h = the angular momentum per unit test mass. For photon, toq = 
=^ /i — >■ 00 and one ends up with the conformally invariant equation but 
without k making its appearance: 



— = -u + 3Mu^ - ^ 
dip-' 2 



= -w + 3Mu^-^. (3) 



In the Schwarzschild-de Sitter (SdS) mcitric. such a disappearance has been 
noted for long [12] but here we find that it occurs despite the presence of 7 
in the metric. The impact of k and 7 on light bending has been investigated 
elsewhere, in Refs. [13-16]. 

Analysis of dynamical system involves converting the second order equation 
into two first order equations [17,18]. For this purpose, we introduce the notation 

d^ fA\ 
u = x, y = x= — (4) 

dip 

to reduce Eq.(2) into a pair of first order autonomous system in the {x,y) phase 
plane 

X = X(x,y) = y (5) 
y = Y{x, y) = a + bx + cx^ + dx~^ + ex~^ (6) 



where 



a=^-lb = -l,c = 3M,d=^,e = -^. (7) 
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First wc discuss stability of circular orbits of light, though their stability is 
not essentially needed. Recall that even in the Schwarzschild spacetime, circular 
light orbits at i? = 3M are unstable. Nonetheless, it is instructive to have a 
look at this aspect in the MKdS solution. 

(a) Massless particle motion 

Light motion occurs in circular orbits defined by R{2 + ^R) — 6M = [see 
Eq.(ll) below] because h? ^ oo, which implies that d = e = but 7^0. 

The equilibrium points are given by i = 0, i/ = 0, which yield ^ ~''^^2c ' 

and ^ -b-s/^^^ iac ^ locate these points on the real phase plane (.t, y), we 

must have = 6^ — 4ac = 1 + Q^M > 0. Now q:^ = => 7 = — g^, so the 
equilibrium points reduce to one single point given by P : (gij,0). For > 0, 
or 7 > — g^, there are two distinct equilibrium points points Q± : (^^^,0) 
where a = + 67M. Thus Q± correspond to two 7— dependent light radii 
i?+ = which expand as follows 

l±Vl+67M ' ^ 



i^. = ^i±^^™^«3M + 0(7), (8) 
7 

i?_=-^-^^™^.-3M-^+0(7). (9) 
7 7 



We have from Eq.(23) below 



<^.=^-'^ (10) 



which yields qo+ = — Vl + 6M7 < leading to unstable radius at i? = R+, 
while go- = \/l + 6M7 > showing that R = R- is a stable radius. The 
basic constraint (reality condition) for both is that 7 > — g^- With, say, 7 = 
—7 X 10~^* cm~^ (value inspired by Ref.[9]), the constraint is always satisfied 
for known lenses (say, for M = 2.9 x lO^^cm, Abell 2744). From Eq.(9), then 
we get the value = 2.86 x 10^^ cm, at which there is stability. However, 
stability of massive particle circular orbits is more relevant, which we examine 
next. 

(h) Massive particle motion 

The equilibrium points are given by i = and y = 0. The equation i; = 
gives r = R = constant, while y = gives 

2MR^ + R\'y-2kR) 

i?(2 + 7i?)-6M • ^ ^ 

The autonomous system (5), (6) can be phrased as a Hamiltonian system as 
follows 

dH 

— = -Y {x, y) = -{a + bx + cx^ + dx~'^ + ex~^) (12) 
BH 

— =X{x,y) = y. (13) 



4 



The necessary and sufRcicnt condition for the system (12), (13) to be a Hamilto- 
nian system, namely, ^ + ^ = 0, is fulfilled for all x and y. Moreover, ^ = 
and therefore H [x, y) = constant (independent of <f). Integrating Eqs.(12),(13), 
we get 

H{x, y) = -{ax + ^x^ + ^x^ - dx'^ - |a;-2) + u{y) (14) 

H{x,y)^]^y^ + v{x) (15) 

where u{y) and v{x) are arbitrary functions subject to the consistency of Eqs.(12) 
and (13). These two equations will match only if 

u{y) = \y^ + C (16) 

v{x) = -{ax + ^x^ + ^x^ - dx-^ - + E (17) 

where C, E are arbitrary constants. The family of Hamiltonian paths on the 
phase plane are given by 

H{x, y) = \y^ - {ax + ^-x^ + ^x^ - dx-' - ^x-^) + G (18) 
where G is a parameter. It follows that 

= -{h + 2cx-2dx-^ -2,ex-^) (19) 
1 (20) 
0. (21) 



dx'^ 



dxdy 

As before, the equilibrium points occur when X = and Y = 0, which give the 
values of the orbit radius r ~ R = constant and as in Eq.(ll). The quantity 
determining stability is [17] 



dx^ dy^ \ dxdy J 



(22) 



Putting the value of h"^, we get at the equilibrium points, suffixed by zero, the 
following expression 

6M R{3kR--,)[R{2 + jR)-6M] 
«° = ^"-^ + R\2kR-j)-2M • ^2^) 

When go > at any point P : {R, 0), we say that P is a stable center, but it 
is an unstable saddle point if g'o < 0. When go = 0, it is an inflection point 
where the system begins to become unstable. Thus go = gives R = -Rmax''') 



5 



beyond which the orbits begin to become unstable. There is also a singular 
radius R = i?sing where go blows up. In all cases we have studied, i?sing < -Re 
for 7 negative and i?sing > -Rds for 7 positive. In either case, the presence of 
singularity is of no concern. 

We shall calculate the values of go for observed lenses. Fortunately, due to 
advances in technology, mass M and Einstein radius Re for several lenses are 
available. The observed lens masses M are used to compute the values of go 
marching from the Einstein radius Re out to Rds- A meaningful stable radius 
should lie between Re and Rds = v^A/S = 1.52 x lO^^cm. Fig.l illustrates for 
one of the observed lenses (Abell 2744, M = 2.90x lO^^cm, Re = 2.97x lO^^cm), 
which shows that stable material radii exist for all radii R > Re when 7 = 
7x 10^^* cm^^ suggesting that the halo might extend even beyond the dS radius, 
which is improbable. Also, the singular radius occurs at i?sing = 8.14 x lO^^cm 
> Rds- (Fig- 2). On the other hand, when 7 — 7 x 10^^* cm^^, stability 
is achieved up to a radius R — 4.25 x 10^^ cm < i?ds beyond which instability 
begins. In this case, the singular radius occurs at i?sing = 9.11 x lO^^cm < Re 
(Fig. 3). Thus the maximum radius up to which stable material circular orbits 
are admissible is i?max'° = 4.25 x 10^^ cm. This limit is not much sensitive either 
to the value of A or to slight magnitude variations in 7. It can be verified that, 
for other lens data (as listed in Ishak et al. [14], but not tabulated here), the 
upper limit remains very nearly the same. 

We conclude the following: A finite upper limit -Rmax^'' the halo radius 
exists within R^s only for a negative 7. We are aware that the tool used here 
is probably far too simplistic than the ground situation, but it is nevertheless 
a good one, which allows us to qualitatively infer a negative 7. Much less is 
yet conclusively known about the halo content or size, which engenders many 
theoretical models, not to mention stiff observational challenges. We considered 
here the MKdS solution by way of illustration but any other metric solution can 
be used likewise. 

One of us (AB) wishes to thank the authorities of the Universita di "La 
Sapienza" , Rome, Italy for financial support during the work. All authors thank 
Guzel N. Kutdusova and Sonali Sarkar for technical assistance. 

Figure captions 

Fig.l. The upper line corresponds to 7 = 7 x 10~^® cm^^ and the lower line 
to 7 = —7 X 10~^* cm~^ in the case of Abell 2744. The plot remains essentially 
the same for other lenses tabulated in Ref.[14]. 

Fig. 2. Curves I correspond to 7 = 7 x 10~^* cm^^ giving a singularity of 
go at i?sing = 8.14 x lO^^cm > R^s- 

Fig. 3. Curves II correspond to 7 = —7 x 10^^^ cm^^ giving a singularity 
of go at i?sing = 9.11 x lO^^cm < Re- 
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Fig.l. The upper line corresponds to 7 = 7 x 10~^^ cm~^ and the lower line 
to 7 = —7 X 10^^^ cm~^ in the case of AbcU 2744. The plot remains essentially 
the same for other lenses tabulated in Ref.[14]. 

Fig. 2. Curves I correspond to 7 = 7 x 10~^^ cm~^ giving a singularity of 
qo at i?sing = 8.14 x lO^^cm > i?dS- 

Fig. 3. Curves II correspond to 7 = — 7 x 10~^* cm~^ giving a singularity 
of qo at iJsing = 9.11 x lO^^cm < R^. 



7 



Figure 1: 



8 



Figure 2: 



9 



Figure 3: 
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